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Abstract 

A new method is suggested to calculate the B — > vrvr hadronic matrix elements 
from QCD light-cone sum rules. To leading order in Us and l/rrib, the sum rule 
reproduces a factorizable matrix element, in accordance with the prediction of the 
QCD factorization approach. Whereas the QCD factorization can only take into ac- 
count the nonfactorizable corrections induced by hard-gluon exchanges, the method 
suggested here also allows a systematic inclusion of soft-gluon effects. In this paper, 
I concentrate on the latter aspect and present a calculation of the nonfactorizable 
soft-gluon exchange contributions to B ^ vrvr. The result, including twist 3 and 
4 terms, is suppressed by one power of l/m-b with respect to the factorizable am- 
plitude. Despite its numerical smallness, the soft effect is at the same level as the 
0(qs) correction to the QCD factorization formula for B vrvr. The method sug- 
gested here can be applied to matrix elements of different topologies and operators 
and to various other i?-decay channels. I also comment on the earlier applications 
of QCD sum rules to nonleptonic decays of heavy mesons. 
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1 Introduction 



Nonleptonic two-body B decays such as B ^ vrvr |l|] are very important tools for studying 
the mechanism of CP-violation (recent overviews of theory and experimental perspectives 
can be found in Refs. 0, 0). In order to fully explore the data on these decays one needs 
reliable theoretical predictions on hadronic matrix elements of the operators entering the 
effective weak Hamiltonian 

Hw = ^VutV:, I (ci(/i) + ^) Oi(/x) + 2c2(/i)Oi(/.) + ...| . (1) 
The operators relevant for B ^ nn are 



0^ = idT^u){uT^b) , Oi = (dT^^u){uT^^b) , (2) 



and the penguin operators denoted by ellipses (for a comprehensive review of the effective 
Hamiltonian see, e.g., Ref. @]). In the above = 7^(1 - 75), Tr(A"A^) = 25"^ and 
O2 = {uT^u){dr'^b) has been Fierz transformed: 

O2 = ^Oi + 2O1 . (3) 

In Eq.(|l|), fi ~ mf, is the normalization scale, and all hard gluon effects at distances 
shorter than are taken into account in perturbative QCD in a form of the Wilson 
coefficients ci,2,...- The remaining contributions at distances > 1/fi have to be included 
in the hadronic matrix elements, such as (vrvr 1 0,1 5). Even an approximate calculation 
of these matrix elements is a very difficult problem for the available QCD methods, in 
particular, for the lattice QCD. Additional complications arise due to the fact that, as a 
rule, several quark topologies (emission, annihilation, penguin, etc.) contribute to each 
matrix element (a general analysis can be found in Ref. [^). 

The method of QCD light-cone sum rules (LCSR) [^, |^, |^ offers a possibility to 
calculate exclusive hadronic amplitudes by matching them, via dispersion relations and 
quark-hadron duality, with the correlation functions of quark currents. These functions 
are evaluated in the spacelike region using OPE near the light-cone. Among other appli- 
cations, various heavy-to-light hadronic form factors have been obtained by this method. 
Importantly, the predicted form factors include both soft (end-point) and hard rescatter- 
ing contributions and reveal a regular 1/mi, expansion (for reviews of the method and 
applications, see Refs. ||, |l^]). In this paper a new LCSR technique is suggested to cal- 



culate hadronic matrix elements for two-body nonleptonic B decays. The sum rule is 
obtained from the light-cone expansion of a three-point correlation function. As a study 
case we consider the B^ —>■ tt+tt" channel. 

One of the main objectives of this work is to quantitatively assess the factorization 
approximation for B vrvr. Recently, the QCD factorization approach has been developed 
|Tl|] for two-body nonleptonic B decays. It was shown that in the limit of the heavy b quark 



mass the decay amplitudes can be treated in QCD using the framework of perturbative 
factorization for exclusive processes [O, with the heavy mass playing the same role as 
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the large momentum transfer. In the hmit mb oo , the effects violating factorization 
in B ^ nn are suppressed either by powers of as{mb) or by powers of KqcD/'i^b- 

The B ^ -R-K amplitude obtained below from LCSR has a very similar structure to the 
one predicted from QCD factorization. It contains a factorizable part, where the B it 
form factor is itself given by LCSR, including both soft and hard contributions. Further- 
more, and it is the main result of this paper, the sum rule approach allows one to calculate 
the effects of soft gluons that violate factorization. In LCSR, these effects are generated 
by higher-twist contributions to the correlation function and are therefore under con- 
trol. Hard nonfactorizable effects correspond to 0{as) corrections to the same correlation 
function, but their estimate demands technically complicated two-loop calculations. The 
new method suggested here avoids certain problems of the earlier applications of QCD 



sum rules to nonleptonic D and B decays |I3|, |I4|, [1^. Moreover, we demonstrate that 
while the use of the light-cone OPE yields a l/rrib behavior of soft higher-twist effects, 
the expansion in local operators employed before fails to reproduce the l/mf, suppression. 

This paper contains a description of the method and the first results for B — >■ tttt. We 
confine ourselves to the contributions of the current- current operators O12 in the emission 
topology. Analysis of other topologies (annihilation, penguin) and operators (quark and 
gluon penguins), as well as applications of the same method to other B decay channels 
will be published elsewhere. 

The content of this paper is as follows. Section 2 contains a derivation of the LCSR for 
B TTTT. In Section 3 the factorizable part of the hadronic matrix element {'n'~7r~^\0i\B'^) 
is reproduced and the soft-gluon nonfactorizable effects determined by (tt^tt+IOiI-B^) 
are calculated. Using these results, in Section 4 we discuss the influence of the soft 
nonfactorizable correction on the QCD factorization formula. Section 5 is devoted to the 
comments on the earlier use of QCD sum rules for nonleptonic decays of heavy mesons. 
A summary is given in Section 6. 



2 The method 

In what follows, we concentrate on the B^ tt+tt^ hadronic matrix elements of the 
operators Oi and Oi. As a starting object for the derivation of LCSR we choose the 
following vacuum-pion correlation function: 

)(p,g,fc) = - | d'xe^^^-^^-J d'ye^^P-'^y{0 \ T{jS(y)O(0)jf ^(a:)} | 7r-(g)) , (4) 

where j^^^ = u'-jal^d and Js^^ = mjji'-y^d are the quark currents interpolating vr and B 
mesons, respectively, and O is either Oi or Oi. The correlator (^ is a function of three 
independent momenta, which are chosen to be g, p — and k. The pion is on-shell, 
= m^. We will work in the chiral limit and set m,r = everywhere, unless there is a 
chirally enhanced combination m^/(m„ -|- md). 

The four-momentum /c 7^ introduced in Eq. @) is unphysical because in reality 
there is no external momentum flow from the weak operator vertex. This momentum 
is added deliberately, otherwise in the correlator the quark states before and after the 
6-quark decay will have one and the same four-momentum, and the dispersion relation in 
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the variable of this momentum squared will contain, below the pole of the ground state 
B meson, a continuum of light "parasitic" contributions. At 7^ 0, the momentum p — q 
in the B channel is independent of P = p — k — q, the total momentum of the state 
formed after the b quark decay. As we shall see below, the procedure is designed in such 
a way that in the final dispersion relation the 4-momentum k automatically vanishes in 
the 5 — i> TTTT ground-state contribution. 

The decomposition of the correlation function in independent momenta is straight- 
forward and contains four invariant amplitudes: 

FP = {P- k)^F(^^ + qjf^ + k^Ff^ + e^,,,q^p'kPFf^ . (5) 

In what follows only the amplitude F^^^ is relevant. The correlation function (^) at A; 7^ 
is, in fact, a 2 — > 2 scattering amplitude. Therefore, counting independent kinematical 
invariants, one has to add two additional variables and p"^ to the four external momenta 
squared = 0, {p — q)^, k"^, and {p — k^. 

The correlation function is calculated in QCD by expanding the T-product of three 
operators, two currents and O, near the light-cone ~ ~ (x — y)^ ~ 0. For this 
expansion to be valid, the kinematical region should be carefully chosen, in order to stay 
far away from hadronic thresholds in both channels of vr and B currents. Consequently, 
the external momenta squared {p — qY and {p — ky, as well as the kinematical invariant P^ 
have to be taken spacelike and large. Simultaneously, to simplify calculations, it is possible 
to put the remaining external momentum squared k'^ and the kinematical invariant p'^ to 
zero having in mind that hadronic states in the channels with four-momenta k and p 
contain heavy b quark and, therefore, the corresponding thresholds are located far enough 
from the zero points. To summarize, our choice of the kinematical region in Eq. is 

q^=p^ = e = 0, \ip-kY\r^\ip-qf\r^\P'\:>AlcD- (6) 

In this region the light-cone expansion of the correlation function (^) is applicable and 
the result is obtained in terms of hard scattering amplitudes convoluted with the pion 
light-cone distribution amplitudes of different twist. The actual calculation is described 
in the next section. At this point, it suffices to represent the result of the QCD calculation 
in a generic form of a dispersion relation in the variable {p — fc)^: 



AO) ,,^2 „2> _ 1 7,. Im^'S'c(s,(p-<^)^^''^ 



Fhc'odP -k)\(p-,)\P') = -j is - • ' . (7) 







To proceed with the derivation, we obtain a corresponding hadronic dispersion relation, 
inserting in Eq. a complete set of intermediate states with the pion quantum numbers: 

(tt) 
h 

where the lowest, one-pion contribution contains the pion decay constant 

{0\u'yal5d\TT{p-k)) = if^{p~k)a, (9) 
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and the hadronic matrix element 



I^^^X{v-q)^P')=^ j rf^xe^(^-«)-(7r-(p-fc)|T{O(0)jf )(x)}|7r-(g)), (10) 

which itself is a two-point correlation function. The spectral function p^^'^^ in Eq. (H) 
accounts for the contributions of excited states and continuum in the pion channel, s^^^^ 
being the lowest threshold. Subtraction terms in both dispersion relations (0) and (P) 
are omitted in anticipation of the Borel transformation which will remove them anyway. 
Physically, the two-point function n^^^ defined in Eq. (|TU]) describes a scattering process. 
In more detail, the pion undergoes a deep elastic scattering due to a combined action 
of the current j'^^ and operator O, with a total spacelike momentum transfer P^. To 
obtain li^^ we follow the usual derivation procedure of QCD sum rules llH] and match 



the dispersion relation (H) with the result (0) of the light-cone expansion at large — A;)^| 
Furthermore, the integral over p|f is approximated using the quark-hadron duality: 

^ P^h\s. {P - q)\ P') _ 1 7. ^^Fqcd{s, {P - g)^ P') 



s — {p — ky IT J s — {p — ky 



where Sq is an effective threshold parameter. After the Borel transformation in the variable 
[p — ky the result reads 

^iZ\iP - q)\ P') = jds e-/^^^Im.F^j^(., {p - q)\ P^) . (12) 



u 

Having at hand an expression for Yi^^J{{p — qY^P'^) valid at large spacelike P^, we 
perform an analytic continuation to large timelike values of P^, keeping the variable {p—qY 
fixed. This procedure is very similar to the evaluation of the timelike asymptotics of the 
pion electromagnetic form factor from a QCD calculation at large spacelike momentum 
transfer. For n^*^^ a natural and convenient point of the analytic continuation is P^ = m\ 
where the invariant mass of two pions is equal to the B meson mass which is sufficiently 
large, ~ m{, ^ ^qcd- The analytic continuation of Eqs. ( |T0| ) and (|12|) yields a relation 



n('^)((p- g)^m|) = zjd'x e'(^-'^)-(7r-(p- fc)7r+(-g)|T{O(0)jf )(x)}|0) 

J ds e-^/'''lm,F^%{s, {p - q)\ m|) (13) 







for the hadronic matrix element which is crossing-related to the one defined in Eq. ([T0|) 
and corresponds to an amplitude of two-pion production with a large invariant mass m\. 
This process is very similar to 7* tt+tt^ but has a more complicated underlying quark- 
current structure. Note also that the analytic continuation to timelike P^ can generate 
in n^'^^((p — g)^, m\) a complex phase, which should be identified, within the accuracy of 
our calculation, with the phase of the final-state rescattering of two pions. 

The next step in the derivation procedure employs analytic properties of the two- 
pion amplitude (^) in the remaining spacelike variable (p — g)^. Inserting in Eq. ([l3|) 
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a complete set of hadronic states with the B meson quantum numbers one obtains the 
following dispersion relation: 



ii,, UP q) , m^) - ^1 _ _ +J s'-{p- qY ' ^ ' 



where the standard definition of the S-meson decay constant {B\j^\()) = m\fB is used. 
Note that in the ground-state contribution the auxiliary momentum k vanishes, due to 
simultaneous conditions {p — q — k)"^ = m\ and {p — g)^ = m^, so that the on-shell 
B ^ Tin matrix element of the operator O is recovered. In the spectral density of excited 
and continuum states k remains nonzero and depends on the value of s' > m^. At large 
— g)^| the relation (|T^) is matched with the result of the QCD calculation given by the 
r.h.s. of Eq. (|13]) which is rewritten in a form of a dispersion relation: 

nf;((p - g)^ = ^J ds e-'/''' J Im.Jm,<J^(., m|).(15) 

b 

The upper limit R of the integration in s' depends, in general, on s and on = m?^. At 
this point one again makes use of quark-hadron duality and approximates the dispersion 
integral over the spectral density p"^'' by the s' > Sq part of the dispersion integral 
(p!5|), where Sq is the effective threshold in the B channel. Finally, the second Borel 
transformation with respect to the variable {p — qY is performed. The resulting LCSR 
for the B^ tt+tt^ matrix element of the operator O reads: 

^ vr+TT-) = {i^-{pW{-q)\0\B{p~q)) 
= , 2 / ^-""'^ I ds' e(-l^-^')/^'^Im., Im.F^J^(., s', m|) , (16) 

U mf 

b 

where R is the upper limit after the duality subtraction in the dispersion integral. 



3 ^ TT+TT hadronic matrix elements 
3.1 Factorizable part 

The remaining task is to calculate Fq^j^ for O = Oi,Oi. At the diagrammatical level, 
there are four topologically different contributions to the correlation function cor- 
responding to four possible combinations of u and d quark-field operators in the pion 
distribution amplitude (0 | Ua{zi)di3{z2) \ n^), where zi = or y, and Z2 = x 01 y 
{a,f3 are the spinor indices). Drawing the quark diagrams we find that each contribu- 
tion yields a i? — *■ tttt matrix element with a certain quark topology: emission {zi = 0, 
Z2 = x), annihilation {zi = 0, Z2 = y), penguin {zi = y, Z2 = x) and penguin annihilation 
[zi = Z2 = y). In what follows we will concentrate on the emission topology adding a 
subscript E to the corresponding quantities. Some of the diagrams contributing to 
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are shown in Figs. 1,2. In the case of the operator Oi there are diagrams, where the 
quarks belonging to the heavy-hght currents do not interact with the quarks of the hght- 
quark currents. The leading-order diagram of this type is shown in Fig. la. To single out 
the contribution of this and all other "factorizable" diagrams, we insert an intermediate 
vacuum state between the weak currents of the operator Oi. Eq. (D is then converted 
into a product of two disconnected two-point correlation functions: 



xlif d^x e^(P-'')^(0 I T{u(0)7^6(0)jf ^(x)} | 7r-(g)) ) , (17) 



where, due to spin-parity conservation, — >• 7^75 (F^ 7^) in the first (second) weak 
current. To proceed further, we insert complete sets of intermediate states with the 
pion and 5- meson quantum numbers in the first and second correlators in Eq. (plT]), 
respectively. After that, the matrix element {'n'~{p)7r^{—q)\0\B{p — q)) is extracted in 
the usual factorized form: 

A^^'\B', TT + TT-) = lUf^M< , (18) 

where /^^ is the B ^ n form factor defined as 

{7r{q)\u^,b\B{p + q)) = 2f+^{p')q, + {f+^{p') + fEAp'))P, ■ (19) 



and taken at = ~ 0. Correction of 0(m^/m^) to Eq. ([18|) is neglected in the chiral 
limit adopted here. 

Within the sum rule method, it is possible to calculate the r.h.s. of Eq. (0) from 
Eq. (p!7D, with a certain accuracy. In leading-order, F^^^ is given by the diagram in 
Fig. la. The answer is easily obtained by forming the pion distribution amplitude 
(0 I Ua{0)dfs{x) I n~), contracting the remaining quark fields in both correlation functions 
and using free-quark propagators. Integrating over the coordinates x,y we obtain: 

,r ^ ■ - i-iP - .)^) / ■ (-) 

In the above, the logarithm corresponds to the bare-loop approximation of the first 
(vacuum-vacuum) correlation function in Eq. (|T^ [|, whereas the integral represents the 
second (vacuum-pion) correlation function. In the latter function, only the leading twist-2 
part is retained, and the pion distribution amplitude (p-j^ is defined in a standard way: 



(0|M(0)7^75rf(x)|7r(g)) = tq,^ rfw e-»«>^(M, /x) , (21) 







where 



(Pn{u, /i) = 6u{l — U) 



n=2,4,. 



(22) 



^The divergent part of the loop vanishes after the Borel transformation and is therefore omitted from 
Eq. (§). 
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p-q k 

(a) 



(b) 



Figure 1: Diagrams corresponding (a) to the leading-order of the correlation function ^ 
for O = Oi; (b) to the higher-twist soft-gluon nonfactorizable contribution for O = Oi. 
Solid, dashed and wavy lines represent quarks, gluons, and external momenta, respec- 
tively. Thick points denote the weak interaction vertices, and ovals the pion distribution 
amplitudes. The cross indicates the point of gluon emission in the second similar diagram. 



C"^"^ are the Gegenbauer polynomials and a„ are multiplicatively renormalizable coeffi- 
cients depending on the normalization scale [i. Following the procedure described in Sect. 
2 and applying, subsequently, duality with Borel transformation in the pion channel, 
analytic continuation — > and duality with Borel transformation in the i?-meson 
channel, one obtains: 







^1,2^7 V' ^.K/^.)J, (23) 



where Uq = ml/sQ and fib is an appropriate scale. In this expression, the first bracket, up 
to a small 0(s/m^) ~ ©(sg/m^) correction, is approximately equal to f-„, because the 
integral over s gives the quark-loop (leading order) term of the SVZ sum rule []TB[ for 



The second bracket in Eq. (p3|) coincides with the leading twist 2 term in the LCSR ^ for 
the form factor f^^. Thus, the factorizable matrix element (|18]) is restored at the leading- 
order level. Factorizable corrections in both correlation functions entering Eq. (O) can 



be added, improving the accuracy of Eq. (|23|). The LCSR for the form factor f^^ [|r 



contains, in addition to the twist 2 term, the twist 3, 4 contributions corresponding to 
two- and three-particle distribution amplitudes of the pion. The numerically important 
twist 3 contribution is proportional to the chirally enhanced coefficient m^/(m„ + md). In 
addition, the 0{as) radiative correction to the twist 2 part has been calculated |jT8|, |T^. 
A certain part of this correction |19| corresponds to the hard rescattering mechanism. 



Summarizing, the hadronic matrix element of the operator Oi in the emission topology 
is given by Eq. ( |T8l) where f^^ itself is calculated from LCSR with an available accuracy. 
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3.2 Nonfactorizable effects 



In the case of the operator Oi, nonfactorizable corrections to the correlation function (^) 
emerge at a two-gluon level and are neglected here. Calculation of the matrix element 
(-7r~7r+|Oi|i?) brings nonfactorizable effects at a one-gluon level into the game. The rele- 
vant correlation function Fjf'^'^ receives contributions of hard gluon exchanges which, at 
0{as), correspond to the diagrams in Fig. 2. These two-loop diagrams can be calculated in 
QCD perturbation theory, but this task is beyond the scope of this paper. In future, it is 
important to obtain this part of the correlation function, not only for the sake of improv- 
ing the accuracy of the LCSR, but also for other important reasons. First, at 0{as), the 
matrix element should partially compensate the scale-dependence of the short-distance 
coefficients ci,2 in the effective weak Hamiltonian. Second, the analytic continuation of 
these two-loop contributions in will generate an imaginary part, yielding a final-state 
rescattering phase in the decay amplitude. Third, it will be interesting to compare the 
0{as) part of the LCSR with the hard nonfactorizable corrections obtained in the QCD 
factorization approach In fact, there is a certain correspondence between the dia- 



grams in Fig. 2 and the quark diagrams contributing to the hard-scattering kernels of 
the QCD factorization. More specifically, the diagrams in Fig. 2a,b, where the hard 
gluon is exchanged between the quarks participating in the weak interaction, are of the 
non-spectator type, whereas the diagrams in Fig. 2c correspond to the hard spectator 
contribution. In particular, it would be instructive to check if at rrih — > oo the form factor 
f^^ factorizes out in the non-spectator part of LCSR, as it does in the QCD factorization 
formula. [| 

There is another class of nonfactorizable effects contributing to F''^^\ with on-shell 
gluons emitted from the quarks of the pion current and absorbed in the pion distribution 
amplitude, as shown in Fig. lb. In terms of the light-cone expansion these contributions 
are of a higher twist, starting from twist 3. The diagrams in Fig. lb are calculated 



employing the light-cone expansion of the quark propagator pO 
^(x,0) ^ -z(0|r{g(x)g(0)}|0) = Jiff 

^ ' 

where G^y = gsG'^^^{\'^ /2), and d is the space-time dimension. In the above, the fixed- 
point gauge for the gluon field is adopted and only the terms proportional to the one 
gluon-field strength are shown, which is the accuracy we need here. 

The relevant quark-antiquark-gluon distribution amplitudes are defined by the matrix 
elements pl | 

{()\u{Q)a^y'-f^Gaf3{vy)d{x)\-K{q)) = ifsnliqaqi^gpu - q^Qt^Qau) 

-{qaqu9p, - qpqu9a,)] J Va, ^3.(a., /i)e-*«(^'"^+^''"«) , (25) 

^Differences may arise due to the fact that in the LCSR approach one of the pions and B meson are 
represented by the integrals over quark spectral densities. In these integrals, quark transverse momenta 
up to 0{sq) and 0{sq ), respectively, are taken into account, whereas in QCD factorization these hadrons 
are described by coUinear light-cone distribution amplitudes. 



9 





Figure 2: Diagrams corresponding to the 0{as) nonfactorizable contributions to the cor- 
relation function ^ for O = Oi. 

{0\u{0)t^,G^p{vy)d{x)\7f{q)) = q^^^^l^^^f^ fva,^\\{a,, fi)e-"^^^^'+y^^'^ 

qx J 

HgiaQp - 9ipqa)f. J Va,^^{a,, ;,)e-^(-^+J/-3)^ (26) 

and the one analogous to Eq. p6| ) with — > 7^,^75, Gap Gap and — > V'lj.x- In the 
above, Gap = ^eappxG^^, Vai = daida2da2,5 (1 - ai - 0^2 - as), and gj^p = gap - {xaqp + 
Xf3qa)/qx. The pion distribution amphtude ^937^ is of twist 3, whereas and ip\\^± are 
of twist 4. In what follows we suppress the n dependence for brevity. 

A straightforward calculation of the two diagrams in Fig. lb yields the following answer 
for the twist 3 contribution: 







{ml - {p- g)2(l - ai)){-P'^va-i - {p - k)'^{l - ^7^3)) 

[(2 -v){q-h)^ 2(1 - v)q ■ {p - k)] {q ■ {p - k)) , (27) 

which is then easily transformed into a form of the dispersion integral (|^): 



jj^iOi) 'rnbf3n f ds f du 



^' ^"'^ 167r2 j s-iv-kf J ml-ip- qfu 
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dv 



V93^(l -U,U-V,V) 



s/{s-P2) 



S + 



mt 



u 



2v 



p2 



:28] 



For simplicity, all terms vanishing after Borel transformations are omitted in this expres- 
sion. In the part of Eq. (pSD which is dual to a pion, the upper limit of the first integral 
is Sq ^ |P^| ~ ml. Hence, the resulting expression for the two-pion matrix element 

Ill^^\{p — q)^,P^) can be simplified by expanding it in powers of s/P^, and neglecting 
very small 0{s/P'^) ~ 0{sq/P'^) corrections: 



n 



(Oi) 

TTTT, tw'd 



{{p-q)\P' 



du 



u{ml — {p — q^u) 



J -u,u-v,v){l + (sl/P') } . 



(29) 



The analytic continuation P^ m\ is then trivial. The twist 4 contribution to n^*^^-* has 
a structure similar to Eq. (pQ]) but a bulky expression not shown here for brevity. Finally, 
the LCSR for the — * tt+tt^ matrix element of the operator Oi is obtained applying to 
n^'^^)((p — g)^, m\) the duality approximation and Borel transformation in the B channel. 
The result is: 



u 

2 

mt 



tm 



B\ 



ds e-'/''" 



( ml 



du 



u 



m\/M''^ -ml/uM'' 



dv 



V?37r(l - U,U - V,V) + f„ 



- 1 



$1(1 



u 



+ 



dv 

V 

s 



( ml 

Um'2 M'2 



3v9j_(l — u,u — v,v) 



1 



1 + , (30) 



(31) 



where the following definitions are introduced: 

^- = ip±{w, 1 - w - v,v) + ip\\{w,l - w - v,v), = $1(1 - v,v) 

The asymptotic form of the pion distribution amplitudes in Eq. (^) is given by 

360aia2a3, (f±{ai) = 10S'^al{l — as), (^||(aj) = — 405^aiQ;2a3 . (32) 



More complete expressions containing scale-dependent nonasymptotic terms can be found 
in Ref. §. 

A few comments are in order concerning LCSR (pO]). Firstly, the magnitude of the 
nonfactorizable effect predicted by this sum rule is determined by the two nonpertur- 
bative parameters /s^r and 6"^. They are defined by the vacuum-pion matrix elements 
{Olua^^j^gsGaisdlTi) and {0\ugsGa^il"'d\7T) , respectively. Secondly, at leading order in 
s^/m^, the matrix element ( pO|) factorizes into and an expression which has a typical 
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structure of the LCSR for a heavy-light form factor. Finally, the analytic continuation 

in does not produce an imaginary part in the nonfactorizable matrix element A^^^\ 
This fact can be interpreted as an absence of a final-state rescattering phase induced by 
soft gluon exchanges. 



3.3 Heavy quark limit 

The heavy-quark mass behavior of the matrix elements (|2^) and (^Oj) can be easily derived. 
To this end, one has to substitute in the LCSR the standard expansions of quantities 
dependent on the heavy quark mass: 

ms = riib + A , Sq =ml + 2mi,UQ , M'^ = 2mfer , = /b/ ^m^, (33) 

where A, ujq, t, and Jb are rrib-independent parameters. One also has to take into account 
the end-point behavior of the light-cone distribution amplitudes, since at — oo, ^ 
1 — ujo/rrif) and only the end-point regions of the momentum fraction u contribute to 
Eqs. (^) and (|0|). As a result we get: 

A^Ub', - vr+vr-) ~ ^(50 ^ vr+vr-) ~ . (34) 

Thus, according to LCSR, the higher-twist nonfactorizable effects are suppressed by one 
power of l/rrif, with respect to the twist 2 factorizable amplitude. 



3.4 Numerical estimates 



The sum rule ( pUD is derived for a finite 6-quark mass, and it is therefore possible to 



obtain a numerical estimate of the nonfactorizable matrix element A^^^\B[ 



and compare it with the factorizable approximation (IT 
To specify the input, we take /^r = 132 MeV, s 



0.7 GeV2 and 



0.5 H- 1.2 



GeV for the parameters of the pion channel. The value of the effective threshold and the 



interval of the Borel parameter are provided by the two-point QCD sum rule |]T6[ for /^r. 
For the parameters of the i?-meson channel and for the pion distribution amplitudes we 
adopt the same input as the one used in LCSR for the form factor f^^ (for a recent update 
see Ref. [||): fe = 180 ± 30 MeV, = 4.7 ^ 0.1 GeV, = 35 ± 2 GeV^ (the values 

obtained from the QCD sum rule for /b), fJ^b = \J^'b ~ '^1 — 2-4 GeV, M'^ 
GeV^, /3^(/ib) = 0.0026 GeV^ 0, H and 5'^{^ib) = 0.17 GeV^ 



10 ± 2 



24, 25 1 . The shapes of 



ISnyH'b) = u.uuzu «orev [J^J, |£0|] auu yiJ^b) 
distribution amplitudes are varied between the BF ansatz of nonasymptotic coefficients 
0] (where, in particular, a2(lGeV) = 0.44, a4(lGeV) = 0.25 and all a„>4 = 0) and purely 



asymptotic form. With this input one, in particular, obtains 



/+^(0) = 0.28 ±0.05. (35) 
To quantify the magnitude of the nonfactorizable soft-gluon effect, we introduce the ratio: 



A^(50^7r+7r-) _4^^)(50 



TUB 



TT + TT ) 



(36) 



12 



and obtain the following estimate 

Ae(52 ^ Tr+vr") = 0.05^0.15 GeV, (37) 

The quoted interval corresponds to the variation of all input parameters within adopted 
ranges and adding linearly the corresponding uncertainties. The parameter A^; is indepen- 
dent of Jb and less sensitive to the input than the individual matrix elements in Eq. (|36D . 
The main uncertainty in A^; is introduced by the nonperturbative parameters /s^ and 5^. 
They are extracted from two-point sum rules with a limited accuracy of about ±30%. 

The soft-gluon contribution to the matrix element (7r"'"7r~|Oi|-B2) turns out to be very 
small, at the level of 1 -i- 2% of the factorizable matrix element (tt+tt^IOiI-B^)- In fact, 
this effect is at the same level as the twist 3,4 contributions of quark- ant iquark-gluon 
distribution amplitudes to the LCSR for the form factor f^^ . 



Soft-gluon correction to the QCD factorization 
formula 



For the 5° tt+tt amplitude, the QCD factorization formula reads ||TT 



X |ci(/x) + + ^C2(/i)F(/z) + O [-^) I , (38) 

where, in order to simplify the discussion, we neglect the matrix elements of penguin 
operators/topologies. The hard nonfactorizable contributions of 0(0;^) parametrized by 



F have been explicitly calculated |]TT|, in a form of the convolution of hard-scattering 
kernels with the pion and 5-meson distribution amplitudes [|. 

The form factor /^^ in Eq. (^) includes not only hard, but also soft contributions 0. 
The latter are not accessible in QCD perturbation theory. Therefore, f^^ represents an 
external input to the QCD factorization formula. Currently, this form factor is calculated 
in lattice QCD using certain extrapolations. It can also be obtained from LCSR, as dis- 



cussed above. The results of both approaches generally agree with each other (see Ref. 
for a recent comparison). The pion distribution amplitudes, at least their nonasymptotic 
terms, e.g., the coefficients a„ in Eq. (0), represent another set of inputs. They, in prin- 



ciple, can be determined from lattice QCD |]5D|. Another possibility to estimate and/or 
constrain the nonasymptotic parts of the pion distribution amplitudes is provided by fit- 
ting the LCSR prediction on the pion electromagnetic form factor to experimental 



data. Yet another external input, used to calculate the hard spectator contribution to F 
in Eq. (^), is the 5-meson distribution amplitude, which is currently taken from models. 

The last term in Eq. ( p8| ) symbolizes contributions of soft nonfactorizable effects. 
Their magnitude cannot be calculated within the QCD factorization approach. The l/m?, 

3 For the B -> TT a similar calculation has been done in Ref. l2ql. 
This important feature distinguishes the QCD factorization approach from other analyses of i? — > tttt 
in perturbative QCD [ETI ESI. 
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suppression of these nonperturbative effects is derived from the infrared behavior of 



perturbative diagrams. Therefore, in order to control the accuracy of QCD factorization at 
finite m^, one has to estimate soft nonfactorizable corrections in a certain nonperturbative 
framework and confirm their l/m^ suppression. 

The LCSR approach described in the previous sections yields for the 
amplitude an expression very similar to Eq. 



A{B', ^ n+n-)LcsR = ^^V^,V:,^[fLmLCSR m||ci(/i) + 



C2(/i) 



+2C2(/i) 



MB', 



'V2 



rriB 



E 



TT ' TT 



rriB 




(39) 



In the nonfactorizable part, in addition to the soft correction in the emission topology 
calculated in Sect. 3, we indicate the contributions to the hadronic matrix element of Oi 
which still have to be calculated. In particular, the terms proportional to Xa,p,pa denote 
soft nonfactorizable corrections due to annihilation, penguin, and penguin annihilation 
topologies. We tacitly assume that all these effects are at least of 0(l/mf,). That has to 
be confirmed by a direct calculation. Finally, 0{as) in Eq. ( p9| ) indicates the contribution 
of hard nonfactorizable effects originating from the diagrams in Fig. 2 and from analogous 
diagrams of other topologies. All these terms can be calculated one by one from the 
underlying correlation function as it was done for . Moreover, the calculation is 
performed using the same framework and input as for the form factor [fB-K{^)\LCSR- We 
stress that in LCSR there is no need to introduce the 5-meson distribution amplitude 
because the heavy meson is interpolated by a current in the correlation function. In 
particular, the hard-gluon exchange with the spectator quark corresponds to the diagrams 
in Fig. 2c, as already explained in Sect. 3.2. 

The completion of Eq. (|39|) is a difficult calculational task. Before it is fulfilled, 
one can use the evaluated ~ XE/rnB term interpreting it as a soft correction to the 
QCD factorization formula, that is, replacing O^AqcD/^b) '^C2Xe/^b in Eq- (pSf). 
To compare the soft and hard nonfactorizable corrections, we take our estimate (|37|) 
and calculate the short- distance part of Eq. (^) numerically. For consistency, we adopt 
fi = fib, adjusting the QCD factorization formula to the characteristic scale at which the 
hadronic matrix elements have been calculated from LCSR. Furthermore, we use ci,2(/^b) 
in NLO and in the NDR scheme, ^^(/if,) = 0.279 corresponding to as{mz) = 0.118, 
and the expression for F given in Ref. [|lT|. In the latter, we adopt for the integral 
Jo ^siO/^ = ^b/Xb over the B meson distribution amplitude the same parameter 

Xb 



0.3 GeV as in Ref. 
in Eq. (|3^ 



11 1. The following estimates are obtained for the separate terms 



Cl(/if, 



C2(/i6 



+ 



a 



— C2{fib)F{^b) 
971 



+ 



2C2(/it 



Xe(B', 



TT ' TT 



niB 



[1.03] + [(-0.007 ^ +0.01) + 0.03i] - [0.005 + 0.015] 



(40) 



where the interval for the real part in the second bracket on r.h.s. corresponds to varying 
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ip.,r{u,Hb) between the BF and asymptotic form in the expression for F. Q This interval 
has only a small correlation with the range in the third bracket on r.h.s. of Eq. (^Ol). The 
latter corresponds to the variation of in Eq. (^) and is to a large extent caused by 
the uncertainty in the nonperturbative parameters /s^ and 6"^. 

We find that the hard and soft nonfactorizable corrections are of the same order 
and their overall effect varies within ±2% of the factorizable amplitude. Thus, in the 

—y TT+TT^ channel both effects are very small and do not influence the decay amplitude, 
which is predominantly factorizable. The situation is different in the case of the C2 + Ci/3 



combination relevant for the colour-suppressed B 



0^0 



TT IT 



decay: 



C2(/ife) + 



Ci(/i6 



+ 



97r 



+ 



2ci(/it 



Xr,(B', 



TT ' TT 



rriB 



[0.103] + [(0.03 ^ -0.04) - 0.104i] + [0.021 ^ 0.064] . 



(41) 



Here the nonfactorizable corrections have a noticeable impact on the amplitude, both on 
its magnitude and phase. 

We conclude that soft nonfactorizable effects are as important as the 0{as) nonfac- 
torizable corrections in the QCD factorization formula. To be cautious, the numbers 
presented in Eqs. (|0]), (|4l|) cannot yet be considered as final predictions to be used in the 
phenomenological analysis. Adding to A^; the remaining soft effects, due to annihilation 
and penguins, will most probably cause further changes in the decay amplitudes. One 
also has to stress that the soft and hard effects in Eqs. and (^) are calculated in 
different approaches. It is more safe and consistent to estimate all nonfactorizable effects 
using one and the same method, and the LCSR approach offers such a possibility. 



5 Earlier applications of QCD sum rules 



There have been several attempts to calculate nonleptonic decays of heavy mesons using 
OPE and QCD sum rules |ll3| , |l^, |15| and focusing on nonfactorizable matrix elements. An 
important prediction of this approach, as emphasized in Ref. |^ , is the channel-dependence 
of nonfactorizable effects. However, certain complications manifested themselves in the 
derivation procedure of the sum rules. The original method of four-point correlation 
functions, applied to two-body D decays, was plagued by a presence of light-hadron states 
in the dispersion relation for the heavy-meson channel. These states, as explained in Sect. 
2, originate due to the absence of an external momentum in the 6-quark decay vertex. A 
different procedure suggested in Ref. [jl4[ for B Dn (and used in Ref. |15| for other chan- 
nels) avoids this problem. If we apply the latter procedure to i? ^ tctt, the starting point is 
then a two-point hadronic matrix element Mbtt = i J d'^x exp(ipx) (7r|T{j^5''(x)O(0)}|-B). 
The expansion of the operator product in a series of local operators. 



(42) 



^We parenthetically notice that the real part of the coefficient F in Eq. (jS^) is quite sensitive to the 
shape of v'7r(M), whereas the imaginary part is independent of this shape, being simply proportional to 
the normalization of the distribution amplitude. 



15 



reduces Mb-k to a set of simpler matrix elements, {i:\Ol! ^\B), which have to be calculated 
separately using, e.g., QCD sum rules. Only the lowest- dimension operator is usually 
taken into account for O = Oi^2- Matching the result of this calculation with the dispersion 
relation for Mbtt in the channel of the current ^^5^ one finally obtains the two-pion 
matrix element (7r7r|0|-B). This estimate is rather ambiguous, because quark-hadron 
duality cannot be applied to the hadronic matrix elements involved in this procedure. 
The method suggested in this paper avoids any unwanted contamination of the dispersion 
relation, due to the auxiliary external four-momentum in the b -decay vertex. With this 
additional element, the sum rule for the hadronic matrix element is directly obtained from 
the initial three-point correlation function, as demonstrated in Sect. 2. One therefore does 
not need to follow the two-stage method of effective operators. 

Apart from these essential improvements, there is also a very important difference 
between LCSR derived in this paper and sum rules obtained before. At all stages of the 
calculational procedure described in Sect. 2,3 we use the light-cone OPE. In contrast, 
the sum rules obtained from four-point correlation functions [1T3 are based on the short- 



distance expansion in local condensates. Similarly, the method of effective operators 
employs a truncated short- distance expansion (|^. This might be a reasonable 



approximation for a specific kinematics of S — > Dtc considered in Ref. where the 



pion energy does not scale with the heavy mass (see also a discussion on this point in 



the second paper of Ref. |TT[). However, as already emphasized in Ref. |jT^, if the non- 
spectator hadron carries a large energy ~ rrib, as in 5 — > vrvr, the resummation of local 
operators with the same twist is necessary. In fact, such a resummation is implicitly 
performed within the procedure of the light-cone OPE employed in Sect. 3. 

It is instructive to demonstrate what happens if one employs a truncated short- distance 
expansion, instead of the light-cone one, in the correlation function (^). We consider 
again the diagram in Fig. lb with the effective operator Oi, and perform a short- distance 



expansion similar to the one in Eq. 

d'y e^(^-'=)^T{jS(2/)Oi(0)} = C^,,,{p - fc)n(0)r^G^^(0)6(0) , (43) 



retaining only one local operator with the lowest dimension. To obtain the corresponding 
Wilson coefficient Ca^iXp one simply has to put t> — > in the argument of the gluon field 
in the quark propagator (^If). This local limit corresponds to the quark propagating 
in the static gluon field, a usual approximation for the condensates in the conventional 



QCD sum rules |]T6|. Thus, in Eq. (^3]) an infinite series of quark-gluon operators of 
the same twist containing derivatives of the gluon field is neglected. Physically, this 
approximation, corresponding to a vanishing four-momentum of the soft gluon, cannot be 
justified, because in the diagram of Fig. lb the momentum of the on-shell gluon is in the 
range < va^q < q, where q is of the order of the virtual quark momenta. 

After substituting the expansion (|43|) back in the correlation function (^, we proceed 
using the light-cone expansion of the product Oeff{0)ji,^\x) and confining ourselves by 
the twist 3 contribution. The result is given by the same expression as in Eq. (pTj), but 
with f = in the denominator: 



short dist. 47i^{-{p - k)^) J ' (mg - (p - g)2(l - a^)) 



16 



X 



3 

-iq-k) + q-{p-k) 



(q-ip-k)). (44) 



We emphasize that in this modified procedure, the short- distance expansion has only been 
partially used. Nevertheless, there are important changes in the resulting expression 



as compared with Eq. (|27|) . The dependence of -F^^^^ on the variable drops out, and 
the dependence on {p — k)"^ reduces to a simple factor l/{p — k^, with an imaginary 
part proportional to 6{s — {p — k)"^). It is therefore difficult to attribute certain duality 
counterparts to the contributions of the pion and hadronic continuum, while matching 

the short-distance version of with the dispersion relation (H). Let us, nevertheless, 
continue the derivation of the sum rule using Eq. (^4]) and assuming a complete pion 
dominance in the dispersion relation. Following the same procedure as described in Sect. 2, 
we obtain the hadronic matrix element 

Ae (5. - " )short dist. = J y ' ' 



X dv V537r(l — U,U — V,V) [ 3 — — f U ) , 



^0 



(45) 



where only the twist 3 part is taken into account. We have checked that a very similar 
result is obtained employing the two-step procedure |]14| based on the effective operators, 
if one uses the approximation ( ^3]) for the expansion (^) and calculates the B ^ it matrix 



element of the effective operator separately, using the LCSR. 

In the heavy quark limit, evaluated using Eq. (|5B|), the r.h.s. of Eq. ( |45| ) is propor- 
tional to m'^/^, enhanced by one power of in comparison with the twist-2 factorizable 
matrix element (PB]). Note that the latter is obtained from the diagram in Fig. la con- 
taining free-quark propagators, and there is no difference between short-distance and 
light-cone expansions in this case. The situation closely resembles QCD sum rules for 
heavy-to-light form factors based on the local condensate expansion. In these sum rules, 
the momenta of soft quarks and gluons forming vacuum condensates are also neglected 
yielding an enhancement of condensate terms at mi, — >■ cxd as compared with the leading- 
order perturbative terms (for more details, see the discussion in Ref. [0]). Numerically, 



the nonfactorizable correction XE/niB calculated from Eq. ( ^5]) is almost an order of mag- 



nitude larger than the estimate (|37D and has a different sign. Note that the sign difference 
in the nonfactorizable matrix element obtained using the new LCSR procedure is a con- 
sequence of the dependence in the correlation function and emerges after the analytic 
continuation (P^ < 0) 



We conclude that even a limited use of the short-distance expansion in the correlation 
function @) yields a sum rule, where it is difficult to separate the ground-state contribution 
applying the usual duality approximation and where the higher-twist terms are enhanced 
at rrib cxD with respect to the leading twist-2 term. 
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6 Conclusions 



In this paper, a new method has been suggested to calculate the hadronic matrix elements 
relevant for two-body nonleptonic B decays. The approach is based on the light-cone ex- 
pansion of a three-point corelation function and on the analytic continuation from the 
spacelike to the timelike region. The resulting hadronic matrix element contains a factor- 
izable part as well as soft and hard nonfactorizable corrections. All these contributions 
can be calculated within one framework using a common input. 

The nonfactorizable soft-gluon exchange correction calculated for the S — > vrvr channel 
is our main result. Importantly, the l/rrif, suppression of this effect has been explicitly 
reproduced, in agreement with the expectation of the QCD factorization approach. More- 
over, it is encouraging that the soft correction in tt^tt~ turns out to be numerically 
small, suggesting that there are no unexpectedly large power corrections to QCD factor- 
ization in this channel. 

The accuracy of the LCSR method is limited, as usual in QCD sum rules, since one 
has to rely on the quark-hadron duality, both in the pion and B channels. Nevertheless, 
one may argue that the duality approximation works reasonably well having in mind, 
e.g. an agreement of the LCSR predictions for the B ^ it form factors with the lattice 
QCD results. It is also important that uncertainties of the sum rule predictions can be 
assessed by varying the input parameters within the currently adopted limits. These 
parameters are universal, therefore, their accuracy can be improved in future by direct 
determinations on the lattice or by fitting QCD sum rules for experimentally measured 
hadronic characteristics to the data. 

The same method, with minor modifications, can be used for all other channels 
where B meson decays either into two light hadrons, or into one charmed hadron (or 
charmonium state) and a light hadron. The important examples in the latter case are 
B — > D^*^7i, J/'ipK^*\ The only limitation is that at least one of the light-quark hadrons 
in the final state has to be a pseudoscalar or vector meson, for which the light-cone 
distribution amplitudes are better known. 
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